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Summary. We investigate the dynamics ol interacting electrons confined to two 
types of quantum dot system, when driven by an external AC field. We first consider 
a system of two electrons confined to a pair of coupled quantum dots jl[ by using an 
effective two-site model of Hubbard-type. Numerically integrating the Schrodinger 
equation in time reveals that for certain values of the strength and frequency of the 
field the tunneling between the dots can be destroyed, thus allowing the correlated 
two-electron states to be manipulated. We then show how Floquet theory |^| can 
be used to predict the field parameters at which this effect occurs. 

We then consider the case of confining the electrons to a single two-dimensional 
quantum dot in the limit of low particle-density. In this system the electrons form 
strongly correlated states termed Wigner molecules, in which the Coulomb inter- 
action causes them to become highly localised in space. Again using an effective 
model of Hubbard-type, we investigate how the AC field can drive the dynamics of 
the Wigner states |H|. As before, we find that the AC field can be used to control the 
tunneling between various charge configurations, and we relate this to the presence 
of avoided crossings in the Floquet quasi-energy spectrum. These results hold out 
the exciting possibility of using AC fields to control the time evolution of entangled 
states in mesoscopic devices, which has great relevance to the rapidly advancing 
field of quantum information processing. 
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1 Introduction 

A quantum dot (QD) is a structure in which electrons can be confined to small 
length scales, comparable to their Fermi wavelength. A set of electrons held in 
such a structure is conceptually similar to a set of atomic electrons bound to 
a nucleus, and for this reason quantum dots are sometimes termed "artificial 
atoms" 0]. Unlike real atoms, the physical properties of quantum dots can 
be easily varied, which gives theorists and experimentalists the opportunity 
to study novel quantum effects in a well-controlled system. 

To extend the atomic analogy further, we can consider linking QDs to- 
gether to form "artificial molecules" . By allowing electrons to tunnel between 
the QDs, the electronic states on the QDs can hybridize, and form new states 
that extend over the whole system || . The degree of the tunneling determines 
the strength of this hybridization. If the tunneling is weak the electrons re- 
main essentially localised on the QDs in analogy with ionic bonding states. 
Conversely, if the tunneling is strong then the electrons form delocalised states 
with a covalent character. Recent transport experiments using AC potentials 
H have been performed on double QD systems, and have indeed revealed 
the ionic or covalent character of the electronic states by measurement of the 
induced photo-current. 

Ever since the pioneering work of Anderson ]7j], it has been known that 
random spatial disorder can cause electronic states to become localised in 
quantum systems. More recently it has been found that an AC driving field 
can produce a similar intriguing effect termed dynamical localisation, in which 
the tunneling dynamics of a particle can be destroyed. One of the first systems 
in which this effect was predicted is that of a particle moving in a double- 
well potential § . A physical realization of this could consist of two coupled 
QDs containing a single electron — the simplest type of artificial molecule 
possible. If this system is prepared with the electron occupying one of the 
QDs, we can expect it to tunnel across to the other QD on a time scale 
set by the Rabi frequency. However, if an AC field of the correct strength 
and frequency is applied to the system, the tunneling is destroyed, and the 
particle will remain trapped in the initial well. 

Weak time-dependent fields are generally treated as small perturbations, 
which produce transitions between the eigenstates of the unperturbed quan- 
tum system. This approach, however, is not applicable to treat the strong 
driving fields required to produce dynamical localisation, and instead the 
technique of Floquet analysis |^] , which is valid in all regimes of driving, has 
proven to be extremely effective. In this approach, which we briefly outline in 
the next section, the important quantities to calculate are the quasi- energies, 
which play a similar role in driven systems to the eigenenergies in the undriven 
case. In particular, dynamical localisation occurs when two quasi-energies of 
states participating in the dynamics approach each other, and become either 
degenerate (a crossing) or close to degenerate (an avoided crossing). Using 
this formalism, analytic and numerical studies of the double- well system have 
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shown P, ^p| , pl| that in the limit of high frequencies, quasi-energy crossings 
occur when the ratio of the field strength to the frequency is a root of the 
Bessel function Jo- 

Adding a second electron to the coupled QD system, however, introduces 
considerable complications. At the low electron densities typically present in 
QDs, strong correlations produced by the Coulomb interaction can signif- 
icantly influence the electronic structure. One of the most dramatic conse- 
quences of this is the formation of Wigner molecule states 1 12 . Understanding 
the interplay between electron correlations and the driving field is, however, 
extremely desirable, as the ability to rapidly control electrons using AC fields 
[ fl3| has immediate applications to quantum metrology H] and quantum 
information processing. In particular, manipulating entangled electrons on 
short timescales is of great importance to the field of quantum computation 

We study this problem here by applying the Floquet formalism to systems 
of interacting particles. The first system we consider is that of two interacting 
electrons confined to a pair of coupled QDs. A consequence of the interaction 
is that the system only responds strongly to the field when the field frequency 
is in resonance with the Coulomb interaction energy. When this condition is 
satisfied we find that, as for the single-particle case, coherent destruction of 
tunneling (CDT) can again occur, but that it governed by the roots of the 
higher-order Bessel functions. We then go on to consider a two-dimensional 
QD in the Wigner regime, which may also be described by a lattice model of 
Hubbard-type Using the same approach we show that CDT can again 
occur when similar conditions are satisfied, and we clarify how an applied 
AC field can drive charge redistributions within a strongly correlated QD. 
Finally we summarize our results and give some brief conclusions. 



2 Methods and approaches 
2.1 Introduction to Floquet theory 

We consider a general quantum system driven by a periodic electric field, 
described by a time-dependent Hamiltonian which we can divide in the fol- 
lowing way: 

H(t)=H t + H I +H AC {t), H AC (t)=H AC (t + nT). (1) 

Here H t holds the tunneling terms, Hj holds the electron-electron inter- 
action terms and H A c{t) describes the interaction of the system with the 
T-periodic driving field. The periodicity of the driving field allows us to 
use the Floquet theorem to write solutions of the Schrodinger equation as 
■0(i) = exp[— iejt](f>j(t) where ej is called the quasi-energy, and <fij(t) is a 
function with the same period as the driving field, called the Floquet state. 
This type of expression is familiar in the context of solid-state physics, where 
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spatial periodicity permits an analogous rewriting of the spatial wavefunction 
in terms of quasi-momenta and Bloch states (Bloch's theorem). 

The Floquet states provide a complete basis, and thus the time-evolution 
of a general state may be written as: 



which is formally analogous to the standard expansion in the eigenvectors of 
a time-independent Hamiltonian. Indeed, in the adiabatic limit, T = 2tt/uj — > 
oo, the quasi- energies evolve to the eigenenergies, and the Floquet states to 
the eigenstates. It is important to note that in this expansion both the basis 
vectors (the Floquet states) and the expansion coefficients explicitly depend 
on time. The nature of this time-dependence is very different however, and 
the superposition of the T-periodicity of the Floquet states with the phase 
factors arising from the quasi-energies produces a highly complicated, quasi- 
periodic time-dependence in general. As the Floquet states have the same 
period as the driving field, they are only able to produce structure in the 
time-dependence on short time-scales. Consequently, the dynamics of the 
system on time-scales much larger than T is essentially determined by just 
the quasi-energies, and hence evaluating the quasi-energies provides a simple 
and direct way of investigating this behavior. In particular, when two quasi- 
energies approach degeneracy the time-scale for tunneling between the states 
diverges, producing the phenomenon of CDT. 

As we shall see for the specific quantum systems we consider, it is fre- 
quently the case that the total Hamiltonian is invariant under the generalized 
parity operation: x — » —x; t — > t + T/2. As a result the Floquet states can also 
be classified into parity classes, depending whether they are odd or even un- 
der this parity operation. Quasi-energies belonging to different parity classes 
may cross as an external parameter (such as the field strength) is varied, but if 
they belong to the same class the von Neumann- Wigner [[l7| theorem forbids 
this, and the closest approaches possible are avoided crossings. Identifying 
the presence of crossings and avoided crossings in the quasi-energy spectrum 
thus provides a necessary (though not sufficient) condition for CDT to occur. 

2.2 Perturbation theory for Floquet states 

Although the quasi-energies are extremely useful for interpretation of the 
time-dependence of a quantum system, they are usually difficult to calculate 
and numerical methods must be employed. When the driving field dominates 
the dynamics, however, it is possible to use a form of perturbation theory 
introduced by Holthaus JTo| | , in which the time-dependent part of the problem 
is solved exactly, and tunneling part of the Hamiltonian, H t , acts as the 
perturbation. This was generalized to treat interacting systems in Refs.[jl]|| 
and was found to be very successful in the high-frequency regime, where hu> 
is the dominant energy-scale. We now give a brief outline of this method. 
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The Floquet states and their quasi-energies may be conveniently obtained 
from the eigenvalue equation: 



where we consider the operator \H(t) — ifid/dt] to operate in an extended 
Hilbert space of T-periodic functions [^8) . The procedure consists of dividing 
the Hamiltonian as in Eq.[l], and finding the eigensystem of the operator 
[Hj + HAc{t) — ifid/dt], while regarding the tunneling Hamiltonian H t as 
acting as a perturbation. Standard Rayleigh-Schrodinger perturbation theory 
can now be used to evaluate the order-by-order corrections to this result, 
requiring only that we define an appropriate inner product for the extended 
Hilbert space: 



Here (-|-) denotes the usual scalar product for the spatial component of the 
wavevectors, and (-|-)t is the integration over the compact time coordinate. 
We shall show in later sections how this method can be used to obtain ana- 
lytical forms which accurately describe the behavior of the quasi-energies for 
the systems we study. 

2.3 Numerical methods 

To study the time-evolution of each system, we used a fourth-order Runge- 
Kutta method to evolve a given initial state in time — typically of the order 
of fifty periods of the driving field. Throughout the time-evolution, physical 
quantities such as the number occupation of a given site were measured, and it 
was ensured that the unitarity of the wavefunction was accurately preserved. 

A number of different methods can be used to numerically calculate the 
quasi-energies of a quantum system, and a detailed description of them is 
given in Rcf. [|| . One technique well-suited to our approach is to evaluate the 
unitary time-evolution operator for one period of the driving field U(t + T,t), 
and then to diagonalize it. It may be easily shown that the eigenvectors of this 
operator are equal to the Floquet states, and its eigenvalues are related to the 
quasi-energies via Xj — exp[— itjT]. This method is particularly convenient 
for our purposes, as C/(T, 0) can be obtained by propagating the unit matrix 
in time over one period of the field, using the same Runge-Kutta method 
described above. 

3 The driven double quantum-dot 




(3) 




We consider a simplified model of this system, in which each QD is replaced 
by a single site. Electrons are able to tunnel between the sites, and we include 
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the effect of interactions by means of a Hubbard-[/ term: 

2 

H = *S ( c i° c 2° + H - c -) + Z (umimi + B *(*H) • ( 5 ) 

0" 2—1 

Here i is the hopping parameter, and for the remainder of this work we 
shall take fi equal to one, and measure all energies in units of t. Ei(t) is 
the external electric potential applied to site i. Clearly only the potential 
difference, E\ — E2, is of physical importance, so we may choose to take the 
symmetric parametrization: 

E E 
Exit) = — coswt, E 2 {t) = —— cos ujt. (6) 

The Hilbert space of Hamiltonian (^) is six-dimensional, comprising three 
singlet states and a three dimensional triplet space. Measurements on semi- 
conductor QDs have shown that the spin-flip relaxation time is typically 
extremely long |l9| ] , and so we have chosen not to include any spin- flip terms 
in the Hamiltonian. Consequently the singlet and triplet sectors are com- 
pletely decoupled, and so if the initial state possesses a definite parity this 
will be retained throughout its time-evolution, and only states of the same 
parity need to be included in the basis. 

To study the time evolution of the system, we used the ground state of the 
static Hamiltonian (a s inglet) as the initial state, and evolved it in time as 



described in Section 22 . Three probability functions were measured through- 
out the time evolution: pLL(t), PRR.(t) and PRL(t), which are respectively the 
probability that both electrons are in the left QD, both are in the right QD, 
and that one electron is in each of the QDs. The Coulomb interaction favors 
separating the electrons, and thus for strong interactions the ground-state 
has a large value of prl, and relatively small values of pll and prr. We 
show in Figjl] the time evolution of these quantities for U = 8 and u> = 4, 
at two different values of electric potential. In both cases the detailed form 
of the time-evolution is highly complicated, but it is clear that the system 
behaves in two distinct ways. In Fig.[j]a the value of prl periodically cycles 
between its initial high value (indicating that each dot holds approximately 
one electron) to nearly zero, while the values of pll and prr correspondingly 
rise and fall at its expense. This behavior is very different to that shown in 
Fig. ^b, where prl never drops below a value of 0.78, and the other two prob- 
abilities oscillate with a very small amplitude. It thus appears that CDT is 
occurring in the second case, and that the system's time evolution is essen- 
tially frozen. We shall term the minimum value of prl attained during the 
time-evolution p m i ni and use this to quantify whether CDT occurs, as a high 
value of p m i n signifies that tunneling has been destroyed, while a low value 
indicates that the electrons are free to move between the QDs. 

In Fig.^Jb we present a contour plot of p m i n as a function of both of the 
frequency and strength of the AC field. Dark areas correspond to low values 
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Fig. 1. Time evolution of the driven double QD system for £7 = 8 and to = 4. (a) 
electric potential, E — 30.0; (b) E = 33.5. Thick solid line = pnL{t), dotted line = 
PLh{t), dashed line = pnii(t). 



of Pmin, and it can be seen that they form horizontal bands, indicating that 
the system is excited strongly by the AC field only at "resonant" values of u>. 
Close examination of this plot reveals that these bands occur at frequencies 
L) = U, U/2, U/3. . ., at which the system can absorb an integer number 
of photons to overcome the Coulomb repulsion between electrons, thereby 
enabling tunneling processes such as | |,J,) — » |0, ||) to occur. We can 
additionally observe that these bands are punctuated by narrow zones in 
which CDT occurs. Their form can be seen more clearly in the cross-section of 
Pmmgiven in Fig.|]a, which reveals them to be narrow peaks. These peaks are 
approximately equally spaced along each resonance, the spacing increasing 
with u>. In Figj^a we show another contour plot of p m i n , this time obtained 
from a full simulation of detailed physical model of two interacting electrons 
confined to a pair of coupled GaAs QDs (for more details on this simulation 
see Ref.|Q]). The striking similarity between these results clearly indicates 
that our simple, effective model (||) indeed captures the essential processes 
occurring in the full system. 

We emphasize that these results are radically different to those obtained 
for non- interacting particles. In this case an analogous plot of delocalisation 
shows a fan-like structure |Tlj ], in which localisation occurs along lines given 
by w = E/xj, where Xj is the j-th root of the Bessel function Jq(x). As 
a test of our method, we repeated our investigation with the inter-electron 
Coulomb repulsion U set to zero, and found that the fan structure was indeed 
reproduced. 

In Fig.||a we show the Floquet quasi-energies as a function of the field 
strength for u = 2, one of the resonant frequencies visible in Fig.||b. We 
see that the system possesses two distinct regimes of behavior, depending on 
whether the driving potential is weaker of stronger than U . For weak fields 
E < U, as studied previously in Ref.j2(|, the Floquet spectrum consists 
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Fig. 2. pmin&s a function of the strength E and energy hco of the AC field: (a) for 
a full simulation of a QD system (huj in units of meV) (b) for the two-site model 
with U — 8 (both axes in units of t). 



of one isolated state (which evolves from the ground state) and two states 
which make a set of exact crossings. Although in this regime p m i„shows lit- 
tle structure, these crossings do in fact influence the system's dynamics. To 
demonstrate this, we show in FigJ^ the Floquet quasi-energies in the weak- 
field regime for the case of U = 16, and plot beneath it the minimum value 
of pll attained during the time-evolution, where this time the state | T4>0) 
has been used as the initial state. It can be seen that for this choice of initial 
condition, the crossings of the quasi-energies again produce CDT and freeze 
the initial state — despite the Coulomb repulsion between the electrons. 

This surprising result may be understood as follows. For large values 
of U , the singlet eigenstates of the undriven system consist of the ground 
state, separated by the Hubbard gap U from two almost degenerate excited 
states. For small values of the driving potential, the two excited states remain 
isolated from the ground state, and constitute an effective two-level system 
with a level-splitting of A ~ 4t 2 /U. Thus if the system is prepared in an 
initial state which projects mainly onto the excited states, its dynamics will 
be governed by the two- level approximation |j 10 llj , and CDT will occur at 



the roots of Jo- We show in Figji^, the quasi-energies obtained from the two- 
level approximation, which give excellent agreement with the actual results 
with no adjustable parameters. As E becomes comparable to the Hubbard 
gap, however, the two excited states are no longer isolated from the ground 
state, and all three levels must be taken into account. This can be seen in the 
progressive deviation of the quasi-energies from the two-level approximation 
as the electric potential approaches U . 

When the electric potential exceeds U, the system displays a very different 
behavior, in which p m i„remains close to zero except at a series of narrow 
peaks, corresponding to the close approaches of two of the quasi-energies. 
A detailed examination of these approaches (see Fig.||b) reveals them to be 
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avoided crossings between the Floquet states which evolve from the ground 
state and the higher excited state, and have the same generalized parity. The 
remaining state, of opposite parity, makes small oscillations around zero, 
but its exact crossings with the other two states do not correlate with any 
structure in p m i n . 




Fig. 3. (a) Quasi-energy spectrum for the two-site model for U = 8 and w = 2, 
circles = exact results, lines = perturbation theory, (b) magnified view of exact 
results for a single avoided crossing. Beneath are the corresponding plots of p m in- 
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Fig. 4. (a) Quasi-energy spectrum for the two-site model for U = 16 and u) = 2, 
circles = exact results, solid line = two-level approximation, e± = ±(A/2)Jo(2E/u>). 
(b) Solid line = minimum value of Pii(t), starting from the initial configuration 
| T-liO). The dotted line denotes p m in, which shows little structure. In both plots 
the vertical dotted line marks the transition to the strong-driving regime. 
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To interpret this behavior in the strong-field regime, we now obtain ana- 
lytic expressions for the quasi-energies via the perturbation theory described 
in Section [2.2| . The first step is to solve the eigenvalue equation (j|) in the 
absence of the tunneling component H t . In a real-space representation the 
interaction terms are diagonal, and so it can be readily shown that an or- 
thonormal set of eigenvectors is given by: 



Imposing T-periodic boundary conditions reveals the corresponding eigenval- 
ues (modulo ui) to be eo = and e± = U. These eigenvalues represent the 
zeroth-order approximation to the Floquet quasi-energies, and for frequencies 
such that U = n u> all three eigenvalues are degenerate. This degeneracy is 
lifted by the perturbation H t , and to first-order, the quasi-energies are ob- 
tained by diagonalizing the perturbing operator Py = {{ei\H t \ej))T- By using 
the well-known identity: 



to rewrite the form of |e±(t)), the matrix elements of P can be obtained 
straightforwardly, and its eigenvalues subsequently found to be eo = and 
e± = ±2J n (E/ui). Fig.||a demonstrates the excellent agreement between this 
result (with n = 4) and the exact quasi-energies for strong and moderate 
fields, which allows the position of the peaks in p m i n to be found by locating 
the roots of J n . Similar excellent agreement occurs at the other resonances. 
For weak fields, however, the interaction terms do not dominate the tunnel- 
ing terms and the perturbation theory breaks down, although we are still 
able to treat the system phenomenologically by using the effective two-level 
approximation. 

4 The square quantum dot 
4.1 What is a Wigner molecule? 

As we remarked in the Introduction, the Coulomb interaction between the 
electrons can significantly affect the electronic structure of a QD. Such strongly 
correlated problems are notoriously difficult to treat, and the addition of a 
time-dependent field complicates the problem even further. When the mean 
inter-electron separation exceeds a certain critical value, however, a surpris- 
ing simplification occurs, as the Coulomb interaction dominates the kinetic 
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energy and drives a transition to a quasi-crystalline arrangement which mini- 
mizes the total electrostatic energy. In analogy to the phenomenon of Wigner 
crystallization in bulk two-dimensional systems (2l]j2^] such a state is termed 
a Wigner molecule |Q. As the electrons in the Wigner state are sharply lo- 
calised in space, the system can be naturally and efficiently discretized by 
placing lattice points just at these spatial locations. A many-particle ba- 
sis can then be constructed by taking Slater determinants of single-particle 
states defined on these lattice sites, from which an effective Hamiltonian of 
Hubbard-type can be generated to describe the low-energy dynamics of the 
system Jl6| . A major advantage of this technique over standard discretization 
[[23|24|] schemes, in which a very large number of lattice points is taken to ap- 
proximate the continuum limit, is that the dimension of the effective Hamilto- 
nian is much smaller (typically by many orders of magnitude) , which permits 
the investigation of systems which would otherwise be prohibitively complex. 
This approach has proven to be extremely successful in treating a variety of 
static problems, including one-dimensional QDs 0, two-dimensional QDs 
with polygonal boundaries ]25| , p6|] , and electrons confined to quantum rings 
[ p7p^| . We further develop this method here by including a time-dependent 
electric field, and study the temporal dynamics of the system as it is driven 
out of equilibrium. 



4.2 Model and Methods 

We consider a system of two electrons confined to a square QD with a hard- 
wall confining potential — a simple representation of a two-dimensional semi- 
conductor QD. Such a system can be produced by gating a two-dimensional 
electron gas confined at a heterojunction interface, and by placing a gate split 
into four quadrants over the heterostructure |||] , the potentials at the corners 
of the QD can be individually regulated. In Fig.[|a we show the ground-state 
charge-density obtained from the exact diagonalization of a square QD |]25f , 
for device parameters placing it deep in the Wigner molecule regime. It can 
be seen that the charge-density is sharply peaked at four points, located close 
to the vertices of the QD. This structure arises from the Coulomb interaction 
between the electrons, which tends to force them apart into diagonally op- 
posite corners of the dot. As there are two such diagonal states, degenerate 
in energy, we can understand the form of the ground-state by considering it 
to be essentially a superposition of these two states (with a small admixture 
of higher energy states). The four points at which the peaks occur define the 
sites on which the effective lattice-Hamiltonian operates, as shown in Fig.||b. 
We take an effective lattice-Hamiltonian of the form: 



H = 



[* (4a c ja + H.c.) + V mnj] n *T™U + E t (t)n t ] . (9) 



Here V represents the Coulomb repulsion between electrons occupying neigh- 
boring sites, and U is the standard Hubbard-t/ term, giving the energy cost 
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(a) (b) 




Fig. 5. (a) Ground-state charge-density for a two-electron square QD. GaAs ma- 
terial parameters are used, and the side-length of the QD is 800 nm, placing it in 
the Wigner regime. The dark areas indicate peaks in the charge-density, (b) Lattice 
points used for the effective lattice-Hamiltonian. 



for double-occupation of a site. As before, Eiit) denotes the electric potential 
at site i, which in general can have a static and a time-dependent component. 
In experiment, static offsets can arise either from deviations of the confin- 
ing potential of the QD from the ideal geometry, or by the application of 
gating voltages to the corners of the QD. Applying corner potentials in this 
way would substantially enhance the stability of the Wigner molecule state, 
and could also be used to ensure that the multiplet of states included in the 
effective lattice-model is well-separated from the other excited states of the 
QD system. In this work, however, we do not explicitly consider the effects 
of static gates, and we neglect the influence of small, accidental offsets en- 
countered in experiment as we expect them to have only minor effects, and 
indeed may even stabilize CDT 30 1. For convenience, we consider applying 
an AC field aligned with the x-axis of the QD, which can be parameterized 
as: 

E E 
Ea = E D = — cosujt, E B — E c = — — cosujt (10) 

where A,B,C,D label the sites as shown in Fig.||b. We emphasize that al- 
though we have the specific system of a semiconductor QD in mind, the 
effective-Hamiltonian we are using can describe a wide range of physical sys- 



tems, including 2x2 arrays of connected QDs |31 , and our results are thus 
of general applicability. 

As for the case of the double QD, we include no spin-flip terms in ([)]) 
and so the singlet and triplet sectors are again decoupled. We choose to 
use initial states with singlet symmetry, which corresponds to the symmetry 
of the system's ground-state. Simple state counting reveals that the singlet 
sector has a dimension of ten, and can be spanned by the six states shown 
schematically in FigJ|, together with the four states in which each site is 
doubly-occupied. 
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Fig. 6. Schematic representation of the two-particle basis states for the singlet sec- 
tor of the Hamiltonian. The ground state of the QD is approximately a superposition 
of states (1) and (2). 

5 Results 



5.1 Interacting electrons, double occupancy excluded 

We begin our investigation by taking the Hubbard-?/ term to be infinitely 
large — that is, to work in the sub-space of states with no double occupation. 
Our Hilbert space is thus six-dimensional, and we use the states shown in 
Fig.^ as a basis. We show in Figj7| the time-dependent number occupation of 
the four sites at two different values of E, in both cases using state (6) as 
the initial state, and setting the AC frequency to u) = 8. In Fig.^a E has a 
value of 100.0, and it can be clearly seen that the electrons perform driven 
Rabi oscillations between the left side of the QD and the right. Accordingly, 
the occupation number of the sites varies continuously between zero and one. 
In Fig.^jb, however, we see that changing the electric potential to a value of 
E = 115.7 produces dramatically different behavior. The occupations of sites 
A and D only vary slightly from unity, while sites B and C remain essentially 
empty throughout the time-evolution. Only a small amount of charge can 
transfer per period of the driving field between the left and right sides of 
the system, producing the small spikes visible in this figure. The amplitude 
of these features is extremely small, however, indicating that the tunneling 
between left and right sides has been almost totally destroyed. 

To confirm that CDT is occurring, we present in Fig]| a comparison of the 
amplitude of the oscillations of ua with the quasi-energy spectrum, as a func- 
tion of the electric potential E. Similarly to the double QD system, we can 
see in Fig.||a that the quasi-energies have two different regimes of behavior. 
The first of these is the weak-field regime, E < V, at which the driving field 
does not dominate the dynamics. In this regime the quasi-energy spectrum, 
and correspondingly, the amplitude of oscillations shows little structure. The 
second regime occurs at strong values of potential, E > V, for which the 
quasi-energy spectrum clearly shows a sequence of close approaches. In Fig.^c 
we show an enlargement of one of these approaches which reveals it to be an 
avoided crossing. Employing the perturbative method described in Section 2.2 
demonstrates that the two quasi-energies involved in these avoided crossings 
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Fig. 7. Time development of the system for U infinite, V — 80 and u> = 8: (a) 
electric potential, E = 100.0 (b) E = 115.7. Solid line indicates the occupation of 
sites A and D, the dotted line the occupation of sites B and C. 



are described by ±2J n (E/uj), where n is equal to V/u>. We may thus again 
think of n as signifying the number of photons the system needs to absorb to 
overcome the Coulomb repulsion between the electrons occupying neighbor- 
ing sites. The results in Fig.|8|b and Fig||d clearly show that the locations of 
the avoided crossings correspond exactly to quenching of the oscillations in 
ha, and so confirm that CDT indeed occurs at these points. 



(a) ' A ' 

/ * A A 






v 

V 








-(b) r 

■ f 

- 

i 

J, 


r 


r 


c 


f]f- 
,1 ■ 


Jf) s 

: \ 


\ 

\j 



50 100 150 200 105 110 115 120 125 
Electric potential, E Electric potential, E 



Fig. 8. (a) Quasi-energies of the system for U infinite, V — 80 and uo = 8: circles 
= exact results, lines=perturbative solution [±2Jio(-E/w)]. (b) Amplitude of oscil- 
lation of the occupation of site A. (c) Detail of quasi-energy spectrum, showing an 
avoided crossing, (d) Detail of amplitude of oscillations. 



5.2 Interacting electrons, double-occupancy permitted 

We now take the most general case, and consider the competition between 
the U and V terms. Setting U to a finite value means that the four doubly- 
occupied basis states are no longer energetically excluded from the dynamics, 
and accordingly we must take the full ten-dimensional basis set. 
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Although it is difficult to obtain precise estimates for the values of pa- 
rameters of the effective Hamiltonian, it is clear that in general U > V. 
Accordingly we choose the parameters U = 160, V = 16 to separate the two 
energy-scales widely for our investigation. We again set the frequency of the 
AC field to uj = 8, and in Figj9|a, we show the quasi-energy spectrum obtained 
by sweeping over the field strength. It is immediately clear from this figure 
that for electric potentials E < U the form of the spectrum is extremely simi- 
lar to the infinite-?/ case. Performing perturbation theory confirms that, as in 
the previous case, the behavior of the quasi-energies is given by ±2J„(E/oj) 
where n = V /uj. We show in Fig.||b the amplitude of the oscillations of tia 
when the system is initialized in state (6), which demonstrates that at the 
locations of the avoided crossings the tunneling parallel to the field is again 
quenched. 
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Fig. 9. (a) Quasi-energies of the system for U — 160 and V = 16, u> = 8: circles = 
exact results, lines=perturbative solution [±2J2(E/ui)]. (b) Amplitude of oscillation 
of the occupation of site A, with (6) as the initial state. 



When the electric potential exceeds the value of U , however, new structure 
appears in the quasi-energy spectrum. A group of four quasi-energies, that 
for weaker fields cluster around zero, become "excited" and make a sequence 
of avoided crossings as the field strength is increased. Perturbation theory 
predicts that these high-field quasi-energies are given by ±2J m (E/cj), where 
m = (U — V)/ui, and thus these avoided crossings arise when the absorption 
of m photons equates to the electrostatic energy difference between the two 
electrons being on neighboring sites, and doubly-occupying one site. This 
then indicates that this structure arises from the coupling of the AC field to 
the doubly-occupied states. 
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To probe this phenomenon, we time-evolve the system from an initial 
state consisting of two electrons occupying site A. In Fig.[iT)|b it can be seen 
that for electric potentials weaker than U the amplitude of the oscillations in 
riA remains small, and shows little dependence on the field. As the potential 
exceeds U, this picture changes, and the AC field drives large oscillations 
in riA, and in fact mainly forces charge to oscillate between sites A and 
B. At the high-field avoided crossings, however, the tunneling between A 
and B is suppressed, which shuts down this process. Instead, the only time- 
evolution that the system can perform consists of undriven Rabi oscillations 
between sites A and D, perpendicular to the field. As these oscillations are 
undriven they have a much longer time-scale than the forced dynamics, and 
thus during the interval over which we evolve the system the occupation of 
A only changes by a small amount, producing the very sharp minima visible 
in Fig. 10 b, centered on the roots of J m {E/u>). 




Fig. 10. (a) Quasi-energies of the system for U = 160 and V = 16, uo — 8: cir- 
cles = exact results, lines=perturbative solution [±2Ji$(E/lo)]. (b) Amplitude of 
oscillation of the occupation of site A, with site A doubly-occupied as the initial 
state. 



As the tunneling perpendicular to the field is undriven, it is straightfor- 
ward to evaluate the time evolution of the initial state, if we assume that 
the left side of the QD is completely decoupled from the right side. The 
occupation of sites A and D is then given by: 

nA(t) — 1 + cos flat, riD — 1 — cos flat (11) 

where J?r = 4i 2 / (U — V) . In Fig.jll] we display the occupations of sites A and 
D as a function of time, for two values of electric potential. At the first value, 
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E = 200, tunneling between the left and right sides of the QD is not quenched, 
and accordingly the occupation of the two sites varies rapidly between zero 
and two as the electrons are driven by the AC field around the system. The 
second value, E = 185.8, corresponds to the first high-field avoided crossing. 
It can be clearly seen that the charge oscillates between sites A and B, with 
a frequency of Hr. These Rabi oscillations are damped, however, indicating 
that the isolation between the left and right sides of the QD is not perfect. In 
this sense we can regard the two sites B and C as providing an environment, 
causing the quantum system composed of sites A and D to slowly decohere 
in time. When the tunneling between the left and right sides of the QD 
is strong, for example at E = 200, this decoherence occurs very rapidly. 
By moving to an avoided crossing, however, and suppressing the tunneling, 
the rate of mixing between the two sides of the QD can be considerably 
reduced, and is just limited by the separation in energy between the two 
quasi-energies. Tuning the parameters of the driving field therefore gives us a 
simple and controllable way to investigate how a two-electron wavefunction 
can decohere in a QD. 




Fig. 11. Time development of the U = 160 and V = 16 system for u> = 8: (a) 
electric potential, E = 200.0 (no CDT) (b) E = 185.8 (CDT). Thick solid line 
indicates the occupation of site A, the thick dotted line the occupation of site D. 
Dotted lines in (b) show the Rabi oscillations of the isolated two-site system, Eq.[ll|. 



6 Conclusions 

In this work we have studied the interplay between Coulomb interactions and 
an AC driving field in two different configurations of QDs. For the case of the 
double QD system we have found that the system presents two different types 
of behavior, depending on whether the driving field dominates the interaction 
energy. For weak driving, we have found the surprising result that a suitable 
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AC field can nonetheless freeze the time-evolution of a doubly-occupied QD 
despite the Coulomb interaction, and that this can be understood by means 
of an effective two-level model. At high field strengths the system is only 
driven strongly by the field at frequencies for which an integer number of 
quanta, n, is equal to the interaction energy. When this condition is satisfied 
CDT is again able to occur at certain well-defined parameters of the field, 
and by using Floquet theory we have shown that these points correspond to 
the roots of J n (E/uj). 

Strong electronic correlations allowed us to use an effective lattice model 
of just four sites to treat the square QD, by taking advantage of the natural 
discretization of the system in a Wigner molecule state. In the effective model, 
the inter-electron Coulomb interaction is described by two parameters, U 
and V, and the dynamics of the system consists essentially of tunneling from 
corner to corner, along the perimeter of the QD. We find again that when 
the frequency of the driving field is in resonance with the Coulomb gap (that 
is, tyllu — V or muj = (U — V)) charge is able to circulate freely around the 
system, except at sharply defined field strengths at which tunneling parallel 
to the field is destroyed. Floquet theory again proved an excellent tool to 
understand this behavior, and revealed that these points correspond to the 
roots of J m {E/uj). 

We have thus shown that AC-fields may not only be used as a spectro- 
scopic tools to probe the electronic structure of QD systems, but can also 
be used to dynamically control the time-evolution of the system. Possible 
applications of CDT range from stabilising the leakage of trapped electrons 
in physical realisations of quantum bits, to acting as "electron tweezers" by 
destroying or restoring tunneling between regions of a mesoscopic device. The 
tunability of the CDT effect, and its ability to discriminate between doubly- 
occupied and singly-occupied states, make it an excellent means for rapid 
manipulation of the dynamics of strongly correlated electrons in mesoscopic 
systems. 
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